Berry phase correction to electron density of states in solids 
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Liouville's theorem on the conservation of phase space volume is violated by Berry phase in the 
semiclassical dynamics of Bloch electrons. This leads to a modification of the phase space density of 
states, whose significance is discussed in a number of examples: field modification of the Fermi-sea 
volume, connection to the anomalous Hall effect, and a general formula for orbital magnetization. 
The effective quantum mechanics of Bloch electrons is also sketched, where the modified density of 
states plays an essential role. 
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Semiclassical dynamics of Bloch electrons in external 
fields has provided a powerful theoretical framework to 
account for various properties of metals, semiconductors 
and insulators I n recent years, it has become increas- 
ingly clear that essential modification of the semiclas- 
sical dynamics is necessary for a proper understanding 
of a number of phenomena. It was known earlier that 
global geometric phase effects 011 Bloch states are 
very important for insulators in our understanding of the 
quantum Hall effect Q, quantized adiabatic purnps [Hi 
and electric polarization |y, L| . It was shown |3, 111 later 
that geometric phase also modifies the local dynamics of 
Bloch electrons and thus affects the transport properties 
of metals and semiconductors. Recently these ideas have 
been successfully applied to the anomalous Hall effect in 
ferromagnetic semiconductors and metals 0,EIH>E1, 
as well as spin transport 0, E3 • 

In this Letter, we reveal a general property of the Berry 
phase modified semiclassical dynamics which has been 
overlooked so far: the violation of Liouville's theorem 
for the conservation of phase space volume. Liouville's 
theorem was originally established for standard classi- 
cal Hamiltonian dynamics, and its importance cannot be 
over emphasized as it serves as a foundation for classical 
statistical physics. The Berry phase makes, in general, 
the equations of motion non-canonical || 0, 0, 0> E3 j 
rendering the violation of Liouville's theorem. Neverthe- 
less, we are able to remedy the situation by modifying 
the density of states in the phase space. 

This modified phase-space density of states enters nat- 
urally in the semiclassical expression for the expectation 
value of physical quantities, and has profound effects on 
equilibrium as well as transport properties. We demon- 
strate this with several examples. First, we consider a 
Fermi sea of electrons in a weak magnetic field, and show 
that the Fermi sea volume can be changed linearly by 
the field. Second, we show how the Berry phase formula 
for the intrinsic anomalous Hall conductivity may be de- 
rived from equilibrium thermodynamics using the Stfeda 
formula [lflj . Third, we provide a general derivation of 
an orbital-magnetization formula which is convenient for 
first-principles calculations. 



k x fl n (k) 
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In addition, we present an effective quantum mechanics 
for Bloch electrons in solids by quantizing the semiclas- 
sical dynamics with the geometric phase. The density of 
states enters in a nontrivial manner into the commuta- 
tors of the phase space coordinates, and relates directly 
to the minimal uncertainty volume in the phase space. 

To begin with, we write down the semiclassical equa- 
tions of motion for a Bloch electron in weak electric and 
magnetic fields |9| 

. _ 1 de n (k) 
r ~ h dk 
hk = -eE(r) -erx B(r) , (lb) 

where Sl n (k) is the Berry curvature of electronic Bloch 
states defined by n„(fe) = i(VfeU n (fe)| X |VfeU n (fe)) with 
\u n (k)) being the periodic part of Bloch waves in the nth 
band, e n {k) is the band energy with a correction due to 
the orbital magnetic moment [see Eq. 1)1 0|l and above]. 
For crystals with broken time-reversal symmetry (such as 
ferromagnetic materials) or spatial inversion symmetry 
(such as GaAs), the Berry curvature fl n (k) is nonzero. 

To show the violation of Liouville's theorem, we con- 
sider the time evolution of a volume element AV = 
ArAfc in the phase space. The equation of motion for 
AV is given by (l/AV)dAV/di = V r • r + V fc • k [IJ. A 
straightforward but somewhat tedious calculation shows 
that the right hand side is equal to —d\n(l+eB-fl/h)/dt, 
which is a total time derivative. Therefore we can solve 
for the time evolution of the volume element and obtain 



AV = AVb/(l + eB ■ £l n /h) 



(2) 



The fact that the Berry curvature is generally k depen- 
dent (and the magnetic field can also depend on r) im- 
plies that the phase space volume element changes during 
time evolution of the state variables (r, k). 

Nevertheless, we have a remedy to this breakdown of 
Liouville's theorem. Equation shows that the volume 
element is a local function of the state variables (through 
the magnetic field and the Berry curvature) and has noth- 
ing to do with the history of time evolution. We can thus 
introduce a modified density of states 



D n (r, k) = (2^)- d (l + eB ■ Sl n /h) 



(3) 
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such that the number of states in the volume element, 
D n (r,k)AV, remains constant in time, where d is the 
spatial dimensionality of the system. The prefactor 
(2-7r)~ d is obtained by demanding that the density of 
states D n (r,k) reduces to the conventional form when 
the Berry curvature vanishes. As will be shown later, 
this density of states corresponds to the minimal quan- 
tum uncertainty volume of the state variables. There- 
fore, it does serve as the semiclassical measure for the 
number of quantum states per unit volume in the phase 
space. Based on this understanding, we write the classi- 
cal phase-space probability density as 



p n (r,k,t) = D n (r,k)f n (r,k,t) , 



(4) 



with f n (r, k, t) being the occupation number of the state 
labeled by (r,fc). Probability conservation demands 
that p n (r,k,t) satisfies the continuity equation in phase 
space. On the other hand, our density of states satis- 
fies dD n /dt = — (V r • r + Vfe • k)D n . It then follows 
that the occupation number introduced above has the 
desired property of being invariant along the trajectory, 
i.e., d/ n /dt = 0. 

We can thus write the real space density of a physical 
observable O in the form [2^ 

0(R)=J2 [ dkD n (r 1 k)f n (r 1 k,t){6S(r-R)) rkn (5) 

where (• • • ) r kn denotes the expectation value in the wave- 
packet state centered at (r, k) with the band index n. In 
the spatially homogeneous case, it reduces to: 



= J2 dkD n (k)f n (k)O n (k) 



(6) 



where O n (k) is the expectation value of O in a Bloch 
state. For simpler notation, we will drop the band index 
n and assume that the integral over k includes the sum 
over n. 

We now discuss the magnitude of the correction term 
eB ■ n/fi to the density of states in Eq. 0. The Berry 
curvature for several materials has been calculated before 
using first-principles method [ill ll2T | . Over large regions 
of the Brillouin zone, its magnitude is on the order of 
a 2 with a being the lattice constant. Thus, eB ~ 
eBa 2 /h is the ratio of the magnetic flux through a unit 
cell to the magnetic flux quantum, and can be 10~ 2 to 
10~ 3 for a magnetic field of 1 tcsla. In the vicinity of 
some isolated points, the Berry curvature can be several 
orders of magnitude higher, leading to bigger effects for 
measurement. In the following, we will present a number 
of applications of our formula Eq(6). 

In our first example, we consider the quantity of elec- 
tron density and show that the Fermi sea volume can 
be changed linearly by a magnetic field when the Berry 



curvature is non-zero. Assuming zero temperature and 
using Eq. ©, we have the electron density as 



dfc 

{2n) d 



1 



eB n 



(7) 



where the upper limit means that the integral is over 
states with energies below the chemical potential p. Not- 
ing that the electron density is fixed by the background 
charge density, we conclude that the Fermi volume must 
change with the magnetic field. To first order, this change 
is given by 



5V F 



^ „ eB n 
dfc 



(8) 



We note that while Landau levels make the Fermi sea 
volume oscillate with the field, the effect described above 
gives an overall shift on average. Such a shift has impor- 
tant implications to those Fermi-surface related behav- 
iors such as transport properties. For instance, in metals, 
it can induce a magnetoresistance linearly depending on 
the magnetic field. On the other hand, in band insula- 
tors, the fc space is limited to the Brillouin zone. Elec- 
trons must populate a higher band if (e/K) J BZ dfc B ■ fl 
is negative. When this quantity is positive, holes must 
appear at the top of the valence bands. Discontinuous be- 
havior of physical properties in a magnetic field is there- 
fore expected for band insulators with non-zero integral 
of the Berry curvatures (Chern numbers). 

In our second example, we show a connection between 
our phase space density of states to the intrinsic anoma- 
lous Hall effect, which is due to spin-orbit coupling in 
the band structure of a ferromagnetic crystal. In the 
context of the quantum Hall effect, Stfeda derived a for- 
mula relating the Hall conductivity to the field derivative 
of the electron density at a fixed chemical potential > 
&xy — — e{dn e / dB z )^. There is a simple justification of 
this relation by a thermodynamic argument by consider- 
ing the following adiabatic process in two dimensions. A 
time dependent magnetic flux generates an electric field 
with an emf around the boundary of some region; and 
the Hall current leads to a net flow of electrons across 
the boundary and thus a change of electron density in- 
side. This argument can be straightforwardly applied to 
the case of anomalous Hall effect and to three dimensions. 
By taking the derivative of the electron density (7J with 
respect to B = Bz at fixed chemical potential, we find 
that 



>xy 



dfc 

(27r) d 



9. z 



(9) 



This is an intrinsic effect because it is independent of 
scattering, and thus differs from conventional skew scat- 
tering and side jump mechanisms 

As a third example of application, we now derive a 
semiclassical formula for orbital magnetization. In the 
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semiclassical picture, a Bloch electron is modeled by a 
wave packet in a Bloch band, which is found to rotate 
about its center of mass in general, yielding an intrinsic 
magnetic moment given by m(k) = — i(e/2?i)(Vfcu]_x 
[Ho(k) — eo(fe)] | Vfcit), where Hq is the Hamiltonian [23j. 
In the presence of a weak magnetic field B, the electron 
band structure energy £o(k) (which may already include 
Zeeman energy from spin magnetization) acquires a cor- 
rection term from this intrinsic orbital moment 0, 0], 
e(k) = £o(fc) — m(k) ■ B. For an equilibrium ensemble of 
electrons, the total orbital magnetization can be found 
from the total energy, which is given by Eq. JBJ as, 



E 



dk 

(27r) d 



1 



iB ft 



eo(k) - m{k) ■ B ) . (10) 



Taking differential of E with respect to B, we obtain the 
magnetization at zero magnetic field to be 



M 



e 

2Ti 



(m(k) -f 

k . j du 
(2Tr) d \dk 



dk 

(2-K) d 
^ dk 



eft 
~~h 

x [2^ 



[Mo - £o(fc)] 

£o(k)-H Q ] 



(11) 

du \ 
dk) ' 



In the upper line of the above expression, the first term 
is the contribution from the intrinsic orbital moment of 
each Bloch electron, and the second term comes from 
the explicit field dependence of the density of states and 
the resulting change in the Fermi volume in Eq. We 
expect this effect to be important in ferromagnetic ma- 
terials with strong spin-orbit coupling. 

Gat and Avron obtained an equivalent result for the 
special case of Hofstadter model [24j . Our derivation pro- 
vides a more general formula that is applicable to other 
systems. Following the discussions on band insulators in 
our first example, there will be a discontinuity of the or- 
bital magnetization if the integral of the Berry curvature 
over the Brillouin zone, or the anomalous Hall conductiv- 
ity, is non-zero and quantized. Depending on the direc- 
tion of the field, the chemical potential /zo in the above 
formula should be taken at the top of the valence bands 
or the bottom of the conduction bands. The size of the 
discontinuity is given by the quantized anomalous Hall 
conductivity times E g /e, where E g is the energy gap. For 
insulators with zero Chern numbers, the orbital magneti- 
zation can be directly evaluated from Wannier functions, 
with result consistent with our general formula [25j . Our 
general formula can also be derived from a full quantum 
mechanical linear response analysis |26|. 

The central result of this paper, equation J3J, can be 
extended to the more general case when Berry curva- 
ture includes the components of ft kr as well as fl kk and 
fl rr 0. In this case, we introduce the Berry curvature 
in phase space, 



ft 



n vr n rk 

ftkr Qkk 



where each block is a 3 x 3 matrix; fl rk = — (fl kr ) T . The 
phase space density of states then reads, 



D = (2tt 



,-d 



det(fi - J) 



(13) 



with J = ( °y J ) . In the special case of electromagnetic 

perturbations with fl kk = e abc fl c , fl^l = -(e/h)e abc B c 
and fl kr = 0, it reduces to J2J). On the other hand, when 
either fl kk or fl rv vanishes, it has a simpler form 



D = {2ir)- d dei(l - Q, 



rk \ 



(14) 



This result has found application in the study of spin- 
force induced charge-Hall effect H?} . 

Finally, we show how the density of states emerges 
naturally in the effective quantum mechanics of Bloch 
electrons. Although our system is not canonical, it can 
nevertheless be quantized following a standard procedure 
developed for non-holonomic systems with second class 
constraints [H HI. First, one redefines the Poisson 
bracket {/,£/}* = {df /d£, a )M ab {dg/d£, b ), where £ a are 
the components of phase space coordinates £ = (r, fe) 
and M = (ft — J) -1 . Our equations of motion can 
then be written as £ a = {£ a ,e}*, where the energy e(£) 
plays the role as the Hamiltonian function. Then, one 
promotes the Poisson brackets into quantum commuta- 
tors: 



{t,£ b } =iM ab 



(15) 



where £ a is the quantum operator corresponding to the 
phase space coordinates. It then follows that a phase 
spa ce p oint acquires a minimal uncertainty volume given 
by " 



i ce p 

m 



= 2 



-d 



det(fi — J) 



-1/2 



(16) 



This can be understood as the phase space volume occu- 
pied by a single quantum state, therefore Eq. H13JI , which 
is proportional to the reciprocal of this volume, can nat- 
urally be regarded as the semiclassical expression for the 
number of quantum states per unit volume in the phase 
space. 

Equation l|15fl presents the effective quantum mechan- 
ics of Bloch electrons. As a demonstration for the validity 
of the quantization scheme as well as the quantum effect 
of the phase space density of states, we consider a simple 
toy model of two dimensional electron system with a con- 
stant Berry curvature, subjected to a uniform magnetic 
field. The commutators read, 



1 



(12) [f >fe ] = K ,^]=i__ 



(17) 
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In the absence of the Berry curvature, we reduce the 
problem to a known case with the familiar nontrivial 
commutator [k x , k y ] = —\(e/Ti)B. In the absence of the B 
field, we have the nontrivial commutator [x, y] — ifi dis- 
cussed extensively in the literature on non-commutative 
geometry. It is interesting to see that in the presence of 
both fields, we do not just have a combination of these 
nontrivial commutators. Instead, we have a nontrivial 
density of states which enters into all of the commuta- 
tors. 

Assuming e(fc) = h 2 k 2 /2m, the system can be solved 
algebraically to yield the energy spectrum and degener- 
acy. We found that the spectrum consists of a set of 
Landau levels with the renormalized cyclotron frequency 
uj c = u}°/[l + (e/U)Bfl\, where w° = eB/m is the usual 



1: 



cyclotron frequency 31]. At the same time, it is more im- 
portant to note that each Landau level still has the same 
degeneracy of eB/h as in the absence of the Berry curva- 
ture. It is known that this degeneracy is directly related 
to the quantized Hall conductance e 2 /h for a filled Lan- 
dau level [13 ■ Had the density of states not entered in 
the commutators, the Landau level degeneracy would be 
modified, violating the topological requirement that the 
Hall conductance for a filled Landau level is quantized. 

Before closing, we note that the phase space density of 
states also enters naturally in the alternative quantiza- 
tion scheme with Feynman path integral. The S matrix 
is calculated by [iil ] 



(out 



S\m) = / IJ[D(€)d£]exp 



Ldt 



(18) 



where L is the Lagrangian for our system || , 



(19) 



with A a {£) = V<>m(0) being the phase space 

gauge potentials associated with the Berry curvature field 

n. 

In summary, we have found a Berry phase correction 
to the phase space density of states for Bloch electrons. 
This correction emerges naturally in both semiclassical 
and quantum mechanics of Bloch electrons, and has pro- 
found effects on the equilibrium and transport proper- 
ties. Because of the fundamental change introduced by 
this correction, it could have important implications on 
other aspects of condensed matter physics, such as the 
Fermi liquid theory. For instance, in the presence of a 
magnetic field, interaction between electrons can change 
the Fermi sea volume by modifying the Berry curvature 
and thus the phase space density of states. 

We acknowledge useful discussions with M. C. Chang, 
R. Resta, G. Y. Guo, Y. G. Yao, and D. Vanderbilt. This 
work is supported by DOE grant DE-FG03-02ER45958. 
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